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Abstract 

Cranked Relativistic Hartree-Bogoliubov theory without and with approx- 
imate particle number projection by means of the Lipkin-Nogami method is 
presented in detail as an extension of Relativistic Mean Field theory with 
pairing correlations to the rotating frame. Pairing correlations are taken into 
account by a finite range two-body force of Gogny type. The applicability 
of this theory to the description of rotating nuclei is studied in detail on the 
example of superdeformed bands in even-even nuclei of the A ~ 190 mass 
region. Different aspects such as the importance of pairing and particle num- 
ber projection, the dependence of the results on the parametrization of the 
RMF Lagrangian and Gogny force etc. are investigated in detail. It is shown 
that without any adjustment of new parameters the best description of exper- 
imental data is obtained by using the well established parameter sets NLl for 
the Lagrangian and DIS for the pairing force. Contrary to previous studies 
at spin zero it is found that the increase of the strength of the Gogny force 
is not necessary in the framework of Relativistic Hartree-Bogoliubov theory 
provided that particle number projection is performed. 
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I. INTRODUCTION 



The development of self-consistent microscopic many-body mean field theories aimed on 
the description of low-energy nuclear phenomena provides necessary theoretical tools for 
an exploration of the nuclear chart into know and unknown regions. This development is 
motivated by theoretical and experimental reasons. Compared with conventional approaches 
such as, for example, the macroscopic-|-microscopic method, self-consistent theories rely on 
a smaller number of assumptions and start from a more microscopic level. For example, 
the starting point of non-relativistic mean field theories is an effective interaction between 
nucleons constituting the nucleus. Such theories based either on zero range Skyrme forces 
or finite range Gogny forces have been widely used starting from seventies. The next step is 
replacing the Schrodinger equation by the Dirac equation and thus considering the relativistic 
mean field (RMF) theory |l|] . In RMF theory, the nucleus is described as a system of point- 
like nucleons, Dirac spinors, which interact in a phenomenological way by the exchange 
of mesons, such as the cr-meson responsible for the large scalar attraction at intermediate 
distances, the cu-meson for the vector repulsion at short distances and the p-meson for the 
asymmetry properties of nuclei with large neutron or proton excess. Such a description has 
a clear advantage that the spin-orbit splitting, which plays an extremely important role 
in low-energy nuclear physics, emerges in a natural way as a genuine relativistic effect. In 
addition, the pseudo-spin symmetry, the origin of which was a long-standing puzzle, does find 
a natural explanation in the framework of RMF theory [Q. RMF theory has been extremely 
successful also in the description of many other facets of low-energy nuclear physics, such 
as ground state properties, giant resonances, superdeformed rotating nuclei etc., see Ref. 
for an overview. 

Since the discovery of the first superdeformed rotational (SD) band in ^^^Dy the 
investigation of superdeformation at high angular momentum remains one of the most chal- 
lenging topics of nuclear structure. The rich variety of physical phenomena at superdeformed 
shapes is based on a complicated and rather subtle interplay of collective and single-particle 
properties. Although at the present stage, a general understanding of this phenomenon has 
been achieved, there are still many unresolved questions related, for example, to the phe- 
nomenon of identical bands and to the treatment of pairing correlations. In addition, the 
microscopic theoretical models used so far are still far from a precise quantitative descrip- 
tion of SD bands which indicates the necessity of further improvements. Different theoretical 
methods based mainly on the concept of the cranking model of Inglis have been employed 
for the quantitative description of various high-spin phenomena at SD shapes, see for exam- 
ple recent overviews in Refs. 0,0. The cranked version of the RMF theory - the Cranked 
Relativistic Mean Field (CRMF) theory [p|- |TT|] is amongst the most successful ones. It has 
been applied in a systematic way to the description of SD bands in different mass regions 
such as A ~ 60 [IIMII, A ~ 80 pl and A ~ 140 - 150 [0,0,|l3-|lll • Pairing correlations 



are expected to be considerably quenched in SD bands of these regions at high spin and thus 
they have been neglected in all studies quoted above. One should clearly recognize that the 
neglect of pairing correlations is an approximation because pairing correlations being weak 
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are still present even at the highest rotational frequencies. Despite this a very successful 
description of many properties of SD bands, such as dynamic J'-^-* and kinematic J^^-* mo- 
ments of inertia, absolute and relative charge quadrupole moments, effective alignments ieff, 
single-particle properties in the SD minimum etc., has been obtained in these studies in an 
unpaired formalism. 

However, the rotational properties of nuclei at low and medium spin are strongly affected 
by pairing correlations. In order to describe such properties within the relativistic framework, 
we have developed the Cranked Relativistic Hartree-Bogoliubov (CRHB) theory. This theory 
is an extension of CRMF theory to the description of pairing correlations in rotating nuclei. 
The brief outline of this theory and its application to the study of several yrast SD bands 
observed in the A ~ 190 mass region has been reported in Ref . |2y] . The present manuscript 
represents an extension of this investigation where both the theoretical formalism and the 
calculations will be presented in much greater details. 

The paper is organized in the following way: In Section |I| a detailed description of 
the CRHB theory without and with approximate particle number projection by means of 
the Lipkin-Nogami method and some of the specific features of the present calculations 
are presented. The shell structure in the A ~ 190 mass region of super deformation and 
the impact of pairing and particle number projection on the rotational and deformation 
properties of rotating nuclei are investigated in detail on the example of the lowest SD 
bands in ^^^Hg and ^^''Pb nuclei in Section |T|. In Section |V|, we study the dependence of 
the results of CRHB calculations with particle number projection on the parametrization of 
the RMF Lagrangian and the Gogny force using SD bands in ^^''Hg and ^^''Pb as an example. 
The properties of yrast SD bands observed so far in even-even nuclei of the A ~ 190 mass 
region are systematically studied in Section 0. Finally, Section Vl summarizes our main 
conclusions. 



II. CRANKED RELATIVISTIC HARTREE-BOGOLIUBOV (CRHB) THEORY 

A. The CRHB equations 

In relativistic mean field (RMF) theory the nucleus is described as a system of point-like 
nucleons, Dirac spinors, coupled to mesons and to the photons. The nucleons interact by 
the exchange of several mesons, namely a scalar meson a and three vector particles u, p and 
the photon. The isoscalar-scalar cx-mesons provide a strong intermediate range attraction 
between the nucleons. For the three vector particles we have to distinguish the time-like 
components and the spatial components. For the photons this means the Coulomb field 
and possible magnetic field in the case where currents play a role. For the isoscalar-vector 
cj-meson the time-like component provides a very strong repulsion at short distances for all 
combinations of particles, pp, nn and pn. For the isovector-vector p-meson the time-like 
components give rise to a short range repulsion for like particles {pp and nn) and a short 
range attraction for unlike particles {up). They also have a strong influence on the symmetry 
energy. In addition, the spatial components of the uo and p-mesons lead to an interaction 
between possible currents, which are for the w-meson attractive for all combinations [pp, 
nn and pn-currents) and for the p-meson attractive for pp and nn-currents but repulsive for 
pra-currents. We have to keep in mind, however, that within mean field theory these currents 
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only occur in cases of time-reversal breaking mean fields as, for instance, in the case of the 
Coriolis fields. 

The starting point of RMF theory is the well known local Lagrangian density 

--n.^n^'' + \mluj,uj^ - \r,.r^'' + Imlp.p^ - \f,^f^'^ (1) 
-gaipaip - gcoi^Y^^t^i^ - gp'^Y^Pt^'ip - eipY^ — A^ijj, 

where the non-linear self-coupling of the cx-field, which is important for an adequate descrip- 
tion of nuclear surface properties and the deformations of finite nuclei, is taken into account 
according to Ref . . The Lagrangian (|l]) contains as parameters the masses of the mesons 



rrio-, ni;^ and nip, the coupling constants g^, g^ and gp and the non-linear terms g2 and g^. 
The field tensors for the vector mesons and the photon field are: 

= d^u^ - d^u^, Rpy = dpPu - dyp^, Fp^ = dpAy - dyAp. (2) 

In the present state of the art of the RMF theory, the meson and photon fields are treated 
as classical fields. 

Two approximations, namely the mean-field approximation, in which the meson field 
operators are replaced by their expectation values 

((t) = (Jo, (^^m) = '^/'Ot^o (3) 

which is justified by the fact that the source terms are large and the No-sea approxima- 
tion, in which only positive-energy states are taken into account , are employed in order 
to solve the Lagrangian (|l]). 

For classical fields we can restrict ourselves to the intrinsic symmetry violating product 
wave function |$) which can be represented as a generalized Slater determinant. As long 
as we consider time-independent static or quasistatic fields, this limits the applicability 
of CRHB theory to the yrast bands and to the excited bands of dominantly quasiparticle 
nature. In particular, the bands based on low-lying collective vibrations, such as for example 
7-, /^-vibrational bands, cannot be described in the present formalism since they have (for 
example, in even-even nuclei) an intrinsic wave function which is a linear superposition 
of many two-quasiparticle states. In order to study such bands one should employ time- 
dependent fields in the random phase approximation in the rotating frame: the task which 
has so far been resolved only in simple non-relativistic models with separable forces 

As discussed in detail in Refs. in the unpaired formalism, the transformation 

to the rotating frame within the framework of the cranking model leads to the cranked 
relativistic mean field equations. Note that in the present investigation we restrict ourselves 
to one-dimensional rotation with rotational frequency VL.j. around the x-axis. Since pairing 
correlations work only between the fermions, the Bogoliubov transformation affects directly 
only the fermionic part. Thus the time-independent inhomogeneous Klein-Gordon equations 
for the mesonic fields obtained by means of variational principle are given in the CRHB 



theory by |20 
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{-A - {n^L^f + ml] a{r) = -g„Ps{r) 



{-A - (l],4)2 + m2}po(r) 
{-A-[fi,(4 + ^,)]' + m2} p(r) 

-A Ao(r) 
-A Afr) 



9pP7{r), 
9pr{r), 

ep'Ar), 
ef{r), 



(4) 



where the source terms are sums of bihnear products of baryon amphtudes 



Ps{r) 


k>0 




Pl'ir) 


fc>0 




pVir) 




- K'{r)M{r), 


rir) 


= 5:[V;"(r)]taV;"(r: 

fc>0 


+ K(r)]tA\4^(r 




k>0 





(5) 



The sums over k > run over all quasiparticle states corresponding to positive energy 
single-particle states {no-sea approximation). In Eqs. (^,^, the indexes n and p indicate 
neutron and proton states, respectively, and the indexes is and iv are used for isoscalar and 
isovector quantities. ffl{r), j^{r) in Eq. (|) correspond to p^^{r) and j^^{r) defined in Eq. 
(I), respectively, but with the sums over neutron states neglected. Note that the Coriolis 
term for the Coulomb potential AQ{r) and the spatial components of the vector potential 
A{r) are neglected in Eqs. since the coupling constant of the electromagnetic interaction 
is small compared with the coupling constants of the meson fields. On this classical level 
only the direct terms of the potentials are taken into account since most parametrizations 
of the RMF Lagrangian have been fitted to experimental data neglecting exchange terms. 
This means that the exchange terms are not fully neglected, but taken into account in an 
averaged way by adjusting the parameters of the direct terms. 

The comparison of Eqs. (QJ^) with Eqs. (4) and (5) in Ref. |TT| indicates that the pairing 
correlations between the fermions have also an impact on mesonic fields through the redef- 
inition of the various nucleonic densities and currents. While the occupation probabilities 
of the single-nucleon orbitals are equal 1 and for the orbitals below and above the Fermi 
level in the systems with no pairing, they are between and 1 when the pairing between 
the fermions is taken into account. 

Contrary to the applications of the Hartree-(Fock)-Bogoliubov theory for the ground 
states of even-even nuclei, the time-reversal symmetry of the intrinsic wave function |$) is 
broken by the Coriolis operator Vt^J^ {Jx is the projection of total angular momentum on 
the rotation axis) in the case of rotating nuclei. Therefore, we do not know a priori the 
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conjugate states in the canonical basis and thus one must solve the full Hartree-Bogoliubov 
problem in the rotating frame |26-28|. The solution of this problem in the RMF theory is to 
large extent similar to the one obtained earlier in the non-relativistic case as far as we treat 
pairing in non-relativistic fashion as it is done in the CRHB theory. Thus in the following 
discussion we will mainly concentrate on the features specific for relativistic case. 

The CRHB equations for the fermions in the rotating frame are given in one- dimensional 
cranking approximation by 





Ek{ ,] (6) 



-A* -h), + \r + nJ: 

where ho is the Dirac Hamiltonian for the nucleon with mass m 

hD = aHV -V{r)) + Voir) + f3{m + S{r)) (7) 

and At- are the chemical potentials defined from the average particle number constraints for 
protons and neutrons (r = p, n) 

{^njN.l^nJ = Z, ($nJiV„|$oJ = N. (8) 

The particle number expectation values {^^iJNt-I^q^) are defined via the (normal) density 
matrices Pp and p„ 

{<l>njNr\<l>nJ=Tr{p,) where p. = ^CK^- (9) 
The Dirac Hamiltonian contains an attractive scalar potential S{r) 

Sir)=gMr), (10) 

a repulsive vector potential Vo(r) 

Voir) = g^uoir) + gpTspoir) + e^-^Ao{r), (11) 
a magnetic potential V{r) 

V{r) = gMr) + 9prMr) + e^^A(r-), (12) 

and a Coriolis term 

— Vt^Jx = ~ ^x{Lx + 2^^)' (''"^) 

The latter two terms are the contributions to the mean field which break time-reversal 
symmetry and induce currents. In the Dirac equation the field V{r) has the structure of a 
magnetic field. Time-reversal symmetry is broken when the orbitals which are time-reversal 
counterparts are not occupied pairwise. At no rotation, this usually takes place in one- 
(multi-)quasiparticle configurations where the magnetic potential V{r) breaks time- reversal 
symmetry. In rotating nuclei the time-reversal symmetry is additionally broken by the 



6 



Coriolis field. The nuclear currents of Eq. (^) are the sources for the space-like components 
of the vector uj{r), p{r) and A{r) fields which give rise to polarization effects in the Dirac 
spinors through the magnetic potential V{r) of Eq. ([I2|) . This effect is commonly referred 
to as nuclear magnetism |^. It turns out that it is very important to take it into account 
for the proper description of currents, magnetic moments p3 and moments of inertia [ITU 



In the present calculations the spatial components of the vector mesons are properly taken 
into account in a fully self-consistent way. 

In Eq. (^), the rotational frequency along the x-axis is defined from the condition |^ 



J=(<l>nJ J. = + (14) 

where / is total nuclear spin. Uk{r) and Vfc(r) are quasiparticle Dirac spinors and E/^ denotes 
the quasiparticle energies. 

The pairing potential (field) A in Eq. (|^) is given by 

A = A„, = lY^VaL^cd (15) 

cd 

where the indices a,b, . . . denote quantum numbers which specify the single-particle states 
with the space coordinates r as well as the Dirac and isospin indices s and r. It contains 
the pairing tensor [] k, 

K = V*U^ (16) 

and the matrix elements V^cd effective interaction in the pp-channel. 

The matrix elements V^^d pairing channel can, in principle, be derived as a one- 

meson exchange interaction by eliminating the mesonic degrees of freedom in the model 
Lagrangian as discussed in detail in the relativistic Hartree-Bogoliubov model developed in 
Ref. |3^. However, the resulting pairing matrix elements obtained with the standard sets of 



the RMF theory are unrealistically large. In nuclear matter the strong repulsion produced by 
the exchange of vector mesons at short distances results in a pairing gap at the Fermi surface 
that is by factor of 3 too large. In addition, standard RMF parameters do not reproduce 
scattering data in the S'o-channel, which is necessary for a reasonable description of pairing 
correlations. On the other hand, since we are using effective forces, there is no fundamental 
reason to have the same interaction both in the particle-hole and particle-particle channel 
| 28| . In a first-order approximation, the effective interaction contained in the mean field F 



is a G matrix, i.e. the sum over all ladder diagrams. On the contrary, the effective force 
in the pp channel (pairing potential A) should be the K matrix, the sum of all diagrams 
irreducible in pp direction. Although encouraging results have been reported in applications 
to nuclear matter [^,^, a microscopic and fully relativistic derivation of the pairing force 



starting from the Lagrangian of quantum hadrodynamics still cannot be applied to realistic 
nuclei. Thus we follow the prescription of Ref. and use a phenomenological interaction 



of Gogny type with finite range (see Eq. ([T7|) below) in the particle-particle channel. Such a 
procedure provides both the automatic cutoff of high-momentum components and, as follows 



^This quantity is sometimes called as abnormal density. 
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from non-relativistic and relativistic studies, a reliable description of pairing properties in 
finite nuclei. 

Although this procedure formally breaks the Lorentz structure of the RMF equations, 
one has to keep in mind that pairing itself is a completely non-relativistic phenomenon. Rel- 
ativistic effects such as the nuclear saturation mechanism due to the cancellation between 
attractive scalar and repulsive vector potentials, spin-orbit splitting and the admixture of 
small components through the kinetic term —icx.'V of the Dirac equation are only im- 
portant for the mean field part of Hartree-Bogoliubov theory and have only a negligible 
counterpart in the pairing field. The pairing density k and the corresponding pairing field 
result from the scattering of pairs in the vicinity of the Fermi surface. The pairing den- 
sity is concentrated in an energy window of a few MeV around the Fermi level, i.e. the 
contributions from the small components of the wave functions to the pairing tensor k are 
very small. Thus in the present version of the CRHB theory, pairing correlations are only 
considered between the positive energy states. Consequently, it is justified to approximate 
the pairing force by the best currently available non-relativistic interaction: the pairing part 
of the Gogny force. At present, this is certainly more realistic than the use of a one-meson 
exchange force in the pairing channel, since this type of interactions has never been opti- 
mized for the description of pairing properties in finite nuclei. Thus the word relativistic in 
CRHB applies only to the Hartree particle-hole channel of this theory. 

The phenomenological Gogny-type finite range interaction is given by 



VPP{1, 2) = f J2 e-^^^'-^''^/"^^' X {Wi + BiP" - HiP^ - MiP^P^) (17) 

=1,2 



where Hi, Wi, Bi, Hi and Mj {i = 1,2) are the parameters of the force and P'^ and P'^ are 
the exchange operators for the spin and isospin variables, respectively. This interaction is 
density-independent. In most of the applications of the CRHB theory in the present article, 
the parameter set DIS |^ (see also Table |) is employed for the Gogny force. Note also that 
an additional factor / affecting the strength of the Gogny force is introduced in Eq. (|T7|). 
This is motivated by the fact that previous studies of different nuclear phenomena within the 
Relativistic Hartree-Bogoliubov theory followed the prescription of Ref. where / = 1.15 



was employed. It turns out, however, that the rotational properties of the SD bands are 
very well described with / = 1.0 if the approximate particle number projection is performed 
by means of the Lipkin-Nogami method. Thus the value / = 1.0 is used in the following if 
no other value is specified. 

In Hartree- (Fock)-Bogoliubov calculations the size of the pairing correlations is usually 
measured in terms of the pairing energy^ defined as 

Epairing = — '^'^T^^^'^^' (^^) 

This is not an experimentally accessible quantity, but it is a measure for the size of the 
pairing correlations in the theoretical calculations. An alternative way to look for the size 
of pairing correlations is to use the BCS-like pairing energies defined as a difference between 



^This quantity is sometimes called as particle-particle correlation energy, see for example Ref. |35] 
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the total energies obtained at a given spin I in the calculations with (£'crhb(-?^)) and without 
{Ecrmf{I)) pairing 

Ebcs{I) — -E'CrhbI-^) — -E'crmf(-^) (19) 

The quantities -Ergs and Epairing coincide only in the BCS approximation. The disadvantage 
of the use of the BCS-like pairing energies is related to the fact that the calculations with 
and without pairing should be performed in order to define this quantity in a proper way. 
The total energy of system in the laboratory frame is given by 

-E'cRHB = E^ + (20) 

where E^ and E^ are the contributions from fermionic and mesonic (bosonic) degrees of 
freedom. Fermionic energies E^ are given by 

E — Epo^j-f + ~l~ Epairing ~l~ -^cm (^1) 

where 

Epart = Tr{hDp) , J = Tr{jxp) , (22) 

are the particle energy and the expectation value of the total angular momentum along the 
rotational axis and 

E^^~fujo^-l^lA-'^^ MeV (23) 

is the correction for the spurious center-of-mass motion approximated by its value in a 
non-relativistic harmonic oscillator potential. 

The bosonic energies E^ in the laboratory frame are given by 

E^^^-^j dr [g, a{r)ps{r) + ^g2a'{r) + Igsa^r)] 
-Ig^ j dr[uJo{r)p^{r)-u,{r)r{r)] 
-\g, j dr[po{v)p\\r) - p{r)r{r)\ 

-le j dr[Ao{r)(f^{r) + A{r)f{r)] 

+ nl jdr [a{r)Llair) - u^,{r)Llu^,{r) + u;(r)(4 + 4)'u;(r) 

-po{r)Llpo{r) + p(r)(4 + S^fp^r)] (24) 

In the systems with broken time-reversal symmetry the spatial components of the cj-mesons 
give larger contributions to the total energy than the ones of the p-mesons because of the 
isovector nature of the p-meson. One should also note that the contribution of the terms 
proportional to to the total energy is very small being typically in the range of ~ 10 — 20 
keV at the highest frequencies of interest in the A ~ 190 mass region and thus, in general, 
can be neglected. 
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B. The Lipkin-Nogami Method 



In recent years it has become clear that a proper treatment of pairing correlations is 
required to describe many nuclear properties. One should note that the Bogoliubov transfor- 
mation is not commutable with the nucleon number operator and consequently the resulting 
wave function does not correspond to a system having a definite number of protons and neu- 
trons. The best way to deal with this problem would be to perform an exact particle number 
projection before the variation However, for heavy nuclei this has been realized so far 



only for non-relativistic separable models [^] in part due to the fact that such calculations 



are expected to be extremely time-consuming for realistic interactions. As a result, an ap- 
proximate particle number projection by means of the Lipkin-Nogami (LN) method [^^-^ 
(further APNP(LN)) is most widely used just because of its simplicity. As illustrated by 
the results of non-relativistic |il-55,^,46,47| and recent relativistic |20,48| calculations, the 



application of this method considerably improves an agreement with experiment especially 
for rotational properties of nuclei compared with the calculations without particle number 
projection. 

Our derivation of the Lipkin-Nogami method is based on the Kamlah expansion |^ 



which is an approximation to the exact projection methods. It is applied to the two-body 
interaction V^'^ in the pairing channel. Note that the Gogny force in the pairing channel is 
density independent in the CRHB theory. For the derivation of the Lipkin-Nogami method 
with density dependent forces see Refs. [^^ . 

In the following discussion, we will limit ourselves to one type of particles. In the numer- 
ical calculations we use this approach for protons and neutrons separately. The eigenstates 
1^"^) of the particle number operator N can be constructed from particle number violating 
Bogoliubov wave functions |$) by means of the particle number projection operator 

1^^) = p^|$) = — d0e*(*-^)'^|$) (25) 
2tc Jo 

where N is actual particle number. With this wave function one obtains the total particle 
number projected energy 



In the case of large particle number and strong deformations in gauge space one expects 
that the Hamiltonian and norm overlap integrals 

/i(0) = (<l>|i^e*^^|<l>) and n{(^) = (<l>|e^^^|<l>), (27) 

are sharply peaked at = and are very small elsewhere. In addition, the quotient 
h{(f)) / n{(f)) should be a rather smooth function. In such a case one can make an expan- 
sion of h{(f)) in terms of n{(j)) in the following way 

M 

h{<j)) = J2 >^mAf"'n{^) (28) 



m=0 



where the Kamlah operator 
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I od) 



($|iV|$) 



(29) 



represents the particle number operator in the space of the gauge angle and are con- 
stants. Mathematically, Eq. ( P^D essentially corresponds to a Taylor expansion of the Fourier 
transformed function /i(0)/n(0) M in Eq. ( pH]) represents the order of expansion, for 

M — >■ cxo this equation is exact. Since we need h{(j)) only in the vicinity of = we can 
already hope to get a very good approximation for this function with a few non-vanishing 
constants Aq, Ai, A2, if they are properly adjusted. The number M certainly depends on 
the widths of the overlap integral which reflects the symmetry violation. The larger 
the symmetry violation, the better the approximation will already be for low M-values. 
Inserting Kamlah operator into Eq. one gets 



M 



h{4>) = ^ A„(<l>|(AAr)V^^|$) with AN = N- (<l>|Ar|$). 



(30) 



m=0 



on Eq. (|8 



The expansion coefficients A^ are determined by applying the operators A/ {k = 0, 1, ...M) 

M 

Sf^h{^) = J2 A™Ar™+'=n(0) , k = O...M. 



m=0 



This yields 



M 



($|(AiV)^i7e^*^|$) = J2 A™($|(AAr)"^+'=e^'^^|$) . 



(31) 



(32) 



m=0 



Taking the limit — this results in a system of the M + 1 linear equations 



M 

A™(AiV)™](AiV)^|<l>) = (33) 

m=0 

for the constants Am- The case of M = 1 corresponds to CRHB theory discussed in Sect. 
11 A| . In the following we will use the shorthand notation (O) = (<l>|0|$). In the case 
of M = 2, the solution of the system of Eqs. (^) under the particle number constraint 
(AN) = leads to the following values of Aq, Ai, A2: 



Ao 
Ai 

Ao 



(H) - \,{{ANf) , 
{HAN) - \2iiANf) 
{(ANy) 

{H[{ANy - {(ANY)]) - {HAN){{ANf)/{{ANY) 



{{AN)^) - ((AiV)2)2 - ((AiV)3)V((AAr)2) 
where the moments of the operator AN are given by 

{{ANf) = 2Tr[x] , 

{{AN)') = ATr[^x] , 

((AiV)^) = ((AiV)2)+8Tr[x(l-6x)], 



(34) 
(35) 

(36) 



(37) 
(38) 
(39) 
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with 



X = p{l-p) and 7 = l-2p. (40) 

Because M is a. hermitian operator with respect to the integration on between and 
27r, one gets a very simple result for the particle number projected energy to order M 



M 

E^ro,= T.^miN-{N)r. (41) 

m=0 



In a full variation after projection method one should vary E^^j.. TheLipkin-No gami method 
consists in treating A2 as a constant during the variation, with its value defined according 
to Eq. (|36D being readjusted self-consistently at each iteration. 
This leads to the variational equation 

^{m- A2(AiV)2|$) - A^($|iV|<|.) = , (42) 

which means that we have to minimize the expectation value of the particle number projected 
energy within the set of product wave functions 1$) under subsidiary condition that A is 
determined by the particle number constraint 

($|iV|$)=Ar, (43) 

provided that the condition A = Ai is fulfilled and the terms proportional to ^A2 are 
neglected. This means in particular that the Lipkin-Nogami method violates the variational 
principle. An extension of the Lipkin-Nogami method to a full variation of the Kamlah 
expansion of the projected energy up to second order is very complicated. So far it has only 



been carried out in non-relativistic theories 50 



The formulation of the Lipkin-Nogami method presented above is given for one kind of 
nucleons (protons or neutrons). In realistic calculations we perform it simultaneously for 
protons and neutrons. It can easily be shown that for M = 2 there is no coupling between 
protons and neutrons. 

The evaluation of the term A2 for a general two-particle interaction V^^ is given as 

1 2TriTri{K,K,~^ V KK~^) — Tr2Tr2{K,*pv an) , , 

4 [Tr(/t/t+)]^ — 2Tr( 

where a = 1 — p and Vabcd = {ab\V^^\cd — dc) is antisymmetrized matrix element of the 
two-particle interaction V^^. The trace Tri represents the summation in the particle-hole 
channel, while Tr2 in particle-particle channel [^]. In Relativistic Hartree-Bogoliubov the- 



ory with the Gogny force in the pairing channel, the sum over the particle-hole part is 
zero, since according to our definition the two-particle interaction V^^ acts only between the 
fermions. In the particle-hole channel we have only classical fields conserving the particle 
number. As a result, the A2 value used in the CRHB calculations with APNP(LN) is given 
by 

_1 Tr2Tr2{K*pvaK) 
4 [Tr(KK+)]2 - 2Tr( 
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Note that in the harmonic oscillator basis used presently in the CRHB(+LN) calculations, 
the density matrix p and pairing tensor k entering into Eq. (|^) are real. 
In Ref. it is shown in detail that Eq. (|45|) can be represented by 



4 [Tr(/t/t+)]^ — 2Tr( 
where the unprojected pairing energy functional is given by 

Ep[K] = ^Tr2Tr2{K*vn). (47) 

Since the modified pairing tensor is much smaller than the pairing tensor k, the main 
contribution to A2 comes from the pairing energy. 

The application of the Lipkin-Nogami method leads to a modification of the Hartree- 
Bogoliubov equations for the fermions, while the mesonic part of the CRHB theory is not 
affected. This modification is obtained by the restricted variation of A2((AA^)^), namely, 
A2 is not varied and its value is calculated self-consistently using Eq. (^61) in each step of 
the iteration. One should note that the form of the CRHB+LN equations is not unique 
(see Refs. PD|J^J^ for details). In the general case, the CRHF+LN equation contains a 
parameter 77 (r^ = 0, ±1) and is given by 




(48) 



where 



hoiv) = hD + 2X2 [(1 + V)P - Trip)] , (49) 

A(r/) = A-2A2(l-r7)/€, (50) 

X{r]) = Ai + A2 [1 + //] , (51) 

Ekiri) = Ek- r/A2 . (52) 

With these definitions and neglecting for simplicity 2A2Tr(p) in ho{'r]) in the following 
discussion it is clear that the case of r/ = +1 corresponds to the shift of whole modification 
into the particle- hole channel of the CRHB+LN theory: hr) h^, + 4A2P leaving pairing 
potential A unchanged. The case of 77 = —1 correspond to the shift of the modification 
into the particle-particle channel A ^ A — leaving ho unchanged. An intermediate 
situation is obtained in the case of r] = 0: ho ^ ho + 2A2P, A — >• A — 2A2fi:. One must note 
that the eigenvalues Ek{rj) of the CRHB+LN equations are identical to the quasiparticle 
energies Ek only in the case of = 0. In the present calculations we are using the case of 
7] = +1 which provides reasonable numerical stability of the CRHB+LN equations. 



C. Physical observables and details of the calculations 

Because, with a few exceptions, the observed SD bands are not linked to the low-spin level 
schemes, the absolute angular momentum quantum numbers / are not known experimentally. 
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Thus the dynamic moment of inertia J^'^^ which contains only the differences A/ = 2 plays 
an important role in our understanding of their structure. In the calculations, the rotational 
frequency VLx^ the kinematic moment of inertia J^^^ and the dynamic moment of inertia j(^) 
are defined by 

Experimental quantities such as a the rotational frequency, the kinematic and dynamic 
moments of inertia are extracted from the observed energies of 7-transitions within a band 



according to the prescription given in Sect. 4.1 of Ref. |5^. One should note that the 
kinematic moment of inertia depends on the absolute values of the spin which in most cases 
are not know in SD bands. 

The charge quadrupole Qq and mass hexadecupole Q40 moments are calculated by using 
the expressions 



IGtt 



go = e^ — V(^'^2o); + 2 {r^Y,,)l (54) 

Q40 = {r^ym)^ + (^'no)„ (55) 

where the labels p and n are used for protons and neutrons, respectively and e is the electrical 
charge. The transition quadrupole moment Qt for a triaxially deformed nucleus is calculated 
by the following expression E 



/I67r / \ cos (7 + 30°) 
V 5 \ 'p cos(30°) 

where the 7-deformation (of the proton subsystem) is defined by 

tan7 = 73TrY (57) 

In the limit 7 — which is typical for SD bands in the A ~ 190 region, Qt is almost 
identical to Qo- The comparison between calculated transition quadrupole moments and 



available experimental data has been performed in Ref. [|20| and thus it will not be repeated 
here. They agree with each other if the uncertainties due to stopping powers are taken into 
account in the experimental data. In addition, it was concluded that more accurate and 
consistent experimental data on Qt is needed in order to carry out such a comparison in 
more detail. 



In the absence of experimentally known spins, an effective alignment approach ||55| , |56|JT8 



plays an extremely important role in the definition of the structure and relative spins of SD 
bands. The effective alignment of two bands (A and B) is simply the difference between 
their spins at constant rotational frequency VLx [^ , ^ , |T^] : 



• B,A 



I. 



(fi,) = /^(fi,)-/^(fi,) (58) 



Experimentally, iejf includes the effects associated with the change of the number of the 
particles and the relevant changes in the alignments of the single-particle orbitals, in the 
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deformation, in pairing etc. between two bands. This physical observable has been used 
frequently in unpaired calculations for the configuration and spin assignments of SD bands 



in the A ~ 60 and A ~ 140 — 150 mass regions (see Refs. |ll2| , p!8| and references therein) 
and for the study of relative properties of smooth terminating bands in the A ~ 110 mass 
region |^. The effective alignment approach exploits the fact that the spin is quantized, 



integer for even nuclei and half-integer for odd nuclei and that it is furthermore constrained 
by signature. One should note that with the configurations and specifically the signatures 
fixed, the relative spins of observed bands can only change in steps of ±2h ■ n, where n is an 
integer number. 

The excitation energies of SD bands relative to the ground state are known definitely 
in ^^"^Hg and ^^^Pb and tentatively in ^^^Pb. The RMF theory with a somewhat simplistic 



pairing and no particle number projection excellently reproduces these data ||58|. It would be 
interesting to see how this result will be changed when the Gogny force is used in the pairing 
channel and APNP(LN) is performed. Considering, however, that such an investigation will 
require a constraint on the quadrupole moment and thus will be extremely time consuming 
within a three-dimensional CRHB computer code, we will leave this question open for future 
studies. 

The CRHB-equations are solved in the basis of an anisotropic three-dimensional har- 
monic oscillator in Cartesian coordinates. The same basis deformation /Sq = 0.5, 7 = 0° and 
oscillator frequency Hujq = 41A~^/^ MeV have been used for all nuclei. All fermionic and 
bosonic states belonging to the shells up to Np = 14 and Nb = 16 are taken into account 
in the diagonalization of the Dirac equation and the matrix inversion of the Klein-Gordon 
equations, respectively. The detailed investigations performed for ^^^Hg indicate that this 
truncation scheme provides reasonable numerical accuracy. The values of the kinematic 
moment of inertia J^^^ and the transition quadrupole moment Qt obtained with such a trun- 
cation of the basis are different from the ones obtained in the fermionic basis with Np = 17 
by less than 1%. The accuracy of the calculated mass hexadecupole moment Q^o is a bit 
lower being ~ 1.5%. The numerical accuracy of the calculations of the physical observables 
within the employed basis has also been tested using different combinations of Po and hujo, 
namely Pq = 0.4, 0.5 and 0.6 as well as Huq = 38^-^/3^ ilA-^Z^ UA-^/\ 47A-^/^, 50^-^/^ 
53 A~^^^ MeV. A similar level of accuracy has been obtained for J^^\ but the new estima- 
tions of accuracy for the calculation of Qt and Q40 are around 2% and 3%, respectively. The 
numerical errors for the total energy (in %) are even smaller. 

At = 0.0 MeV, the single-particle orbitals are labeled by means of the asymptotic 
quantum numbers [A^n^Ajfi (Nilsson quantum numbers) of the dominant component of the 
wave function and superscripts to the orbital labels (e.g. [651]l/2+) are used to indicate the 
sign of the signature r for that orbital (r = 



III. THE IMPACT OF PAIRING AND PARTICLE NUMBER PROJECTION 

In the present section, the impact of pairing and approximate particle number projection 
will be studied in detail on the example of the yrast SD bands observed in ^^^Hg and 
^^"^Pb. In order to differentiate different types of calculations, the following abbreviations 
are introduced: 

• CRMF - cranked relativistic mean field calculations without pairing correlations. 
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• CRHB - cranked relativistic Hartree+Bogoliubov calculations without particle number 
projection, 

• CRHB+LN - cranked relativistic Hartree+Bogoliubov calculations with approximate 
particle number projection by means of the Lipkin-Nogami method (APNP(LN)). 

In the following the lines representing the results of different CRHB(+LN) calculations will 
be labelled in the figures by 

RMFset + / * Gogset + LN 

where RMFset and Gogset are the RMF and Gogny forces used in the calculations, / 
is the scaling factor for the Gogny force (see Eq. ([T7|) ). Note that / will in general be 
omitted when / = 1.0. LN is used only in the cases when APNP(LN) is performed in the 
calculations. The lines showing the results of the CRMF calculations will be labelled by 
RMFset [CRMF]. 

We start from the results of the CRMF calculations. Neutron and proton single-routhian 
diagrams calculated with the parameter set NLl (see Table H) and drawn along the defor- 
mation path of the yrast SD configurations in these two nuclei are shown in Fig. |l|. The 
shell structure of ^^^Hg is dominated by the large Z = 80 and = 112 SD shell gaps. The 
same shell gaps appear also in the calculations with the set NL3, see Fig. 0, but compared 
with the set NLl the Z = 80 SD shell gap is somewhat smaller while the N = 112 SD shell 
gap is more pronounced. Similar shell gaps are seen also in ^^^Pb (see right panels of Fig. 
|I]). Considering that the Z = 82 SD shell gap decreases with increasing rotational frequency 
due to down-sloping 7r[651]l/2+ orbital, it is more reasonable to consider the nucleus ^^^Hg 
as a doubly magic in this mass region. While the NLl and NL3 forces give the same set of 
single-particle orbitals in the vicinity of magic SD shell gaps, their relative ordering (ener- 
gies) is somewhat different. As studied in detail in Ref. |18[ (Section 6.4), the difference in 
the single-particle energies at superdeformation when different RMF parametrizations are 
used is to a great extend connected with their energy differences at spherical shape. Note 
that the proton single-particle spectra in the vicinity of the Z = 80 SD shell gap in the 
A ~ 190 mass region reveal large similarities with the neutron single-particle spectra in the 
vicinity of the A^ = 80 SD shell gap in the A ~ 140 — 150 mass region, see Fig. ^ in the 



present article and Fig. 17 in Ref. ||18|| . 



Figs. |l| and show some similarities and differences with the results obtained in the 
non- relativistic calculations with the Wood-Saxon potential (see Fig. 3 in Ref. [^), Skyrme 
forces (see Figs. 6 and 7 in Ref. |4l|]) and Gogny forces (see Fig. 5 in Ref. [Q). Both 



in relativistic and non-relativistic calculations a similar shell structure and a similar set of 
single-particle states (although their relative energies are different) appear in the vicinity of 
the SD shell gaps. The main difference with non-relativistic calculations is related to the 
larger size of the SD shell gaps and to the lower level density in the vicinity of these gaps. 
These features are connected with low effective mass of RMF theory, see the discussion in 



Ref. for details. 

As can be seen in Fig. |^, the CRMF calculations do not reproduce the experimental 
kinematic and dynamic moments of inertia. In the case of ^^^Hg, the calculated values of J'-^-' 
and J*^^-* are almost constant, while in the case of ^^^Pb the unpaired proton band crossing is 
clearly seen. It originates from the interaction between 7r[642]5/2+ and 7r[651]l/2+ orbitals. 
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see bottom panels in Fig. |l|. The orbital 7r[642]5/2"'" is occupied before band crossing in 
the lowest SD configuration, while the orbital vr [651] 1/2+ is occupied after band crossing. 
Note that in the A ~ 150 mass region, the crossings observed in some bands of the nuclei 
around ^^^Gd have been attributed to the unpaired interaction of these two orbitals, but in 
the neutron subsystem, see Refs. [|rT|,^. 

The inclusion of pairing without particle number projection somewhat improves the 
agreement with experiment, see Fig. ^. Indeed, the slope of the calculated kinematic and 
dynamic moments of inertia at low rotational frequencies is coming closer to experiment, 
but still the disagreement is considerable. In addition, a proton pairing collapse takes place 
in the CRHB calculations. This collapse is calculated in ^^^Pb at Qx ~ 0.35 MeV and 
in ^^^Hg at ~ 0.10 MeV, see Figs. ^ and In the case of ^^^Pb it is correlated 
with the alignment of the 7r[651]l/2"'" orbital which reveals itself in a sharp increase of 
J(2), see Fig. |a. This crossing coincides with the one seen in the CRMF calculations. 
At higher frequencies, the calculated kinematic and dynamic moments of inertia for the 
proton subsystem and the transition quadrupole moments obtained in the CRMF and CRHB 
calculations are very close, see Figs. |^ and ^,d, respectively. However, at lower frequencies 
proton pairing is still important and the results of the calculations with and without pairing 
for proton J^^^ and J*-^-* values are different. Contrary to the proton subsystem, the pairing 
correlations in the neutron subsystem do not collapse in the rotational frequency range 
under consideration. Neutron pairing energies E^^-^-^g, being larger in absolute value than 
the proton ones at fl^ = 0.0 MeV, start around —6 MeV and then smoothly decrease in 
absolute value with increasing rotational frequency coming close to MeV at Qx = 0.5 
MeV, see Figs. ^,b. At low and medium rotational frequencies, there is a considerable 
difference between the corresponding neutron moments of inertia (J*-^-* or J*^^-*) obtained 
in the CRMF and CRHB calculations, see Fig. ^. This difference, however, is small at 
the highest rotational frequencies, where the magnitude of neutron pairing correlations is 
negligible. 

Note that neutron pairing correlations alone do not lead to a considerable modification of 
transition quadrupole moments Qt compared with the CRMF calculations, see an example 
of ^^^Hg (Figs. P3,d). On the contrary, in ^^^Pb, where both neutron and proton pairing 
persist up to Qx ~ 0.35 MeV, the transition quadrupole moment Qt obtained in the CRHB 
calculations is larger than the one of the CRMF calculations. 

The inclusion of approximate particle number projection by means of the Lipkin-Nogami 
method leads to a very good agreement between the calculated and experimental kinematic 
and dynamic moments of inertia, see Fig. ^. The correlations induced by the Lipkin-Nogami 
prescription compared with CRHB produce the desired effects, namely, to diminish the kine- 
matic moments of inertia at all frequencies, to diminish the dynamic moments of inertia at 
low and medium rotational frequencies and to delay proton and neutron alignments to higher 
frequencies. These effects are due to stronger pairing correlations seen in CRHB-I-LN com- 
pared with CRHB (Figs. §a,b). Contrary to the CRHB approach where the proton pairing 
collapses at some frequencies, no such collapse appears in the CRHB+LN calculations in the 
rotational frequency range under consideration. Thus the proton and neutron subsystems re- 
main correlated even at Qx = 0.5 MeV. The decrease of the magnitude of pairing correlations 
with increasing rotational frequency seen both in CRHB and CRHB+LN is due to a strong 
Coriolis anti-pairing effect. The transition quadrupole moments Qt calculated in CRHB+LN 
show a different behavior as a function of rotational frequency compared with CRMF and 
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CRHB, see Figs. ^,d. While in the latter approaches, the Qt values decrease with increasing 
rotational frequency, they initially increase and then decrease in the CRHB+LN approach. 
The calculated decrease of Qt is due to an anti-stretching effect caused by the Coriolis force, 
but the differences in the behavior of Qt as a function of VLx should be attributed to different 
strengths of the pairing correlations in the CRHB and CRHB+LN calculations. 

The increase of the kinematic and dynamic moments of inertia is a complex effect which 
predominantly includes the gradual alignments of the pairs of ji5/2 neutrons and ^13/2 protons 
and the decreasing pairing correlations with increasing rotational frequency. The quasipar- 
ticle routhian diagrams of Fig. |] shows the typical quasiparticle spectra obtained in the 
CRHB+LN calculations and the alignments in the above discussed pairs. Note that the 
alignment of the proton pair takes place at a higher rotational frequency than the one of the 
neutron pair. It is reasonable to expect that the behavior of the total dynamic moment of 
inertia will sensitively depend on the balance of the alignments in the proton and neutron 
subsystems. Indeed, the total J*^^) in ^^^Pb does not show the decrease above the frequency 
of full alignment of the neutron pair, while such a decrease is seen in ^^^Hg (see Figs. |^,c). 

Total (proton + neutron) BCS-like pairing energies -Ebcs (see Eq. (|T^)) obtained for the 
lowest SD configurations in ^^^Hg and ^^^Pb in the calculations with and without APNP(LN) 
are shown in Fig. p. These quantities have more physical content than Epairing defined in 
Eq. (0) since they directly show the gain in binding energy due to the pairing correlations. 
The comparison with Figs. ^,b allows to conclude that -Ergs is smaller by roughly an 
order of magnitude than Ep°J;-^-^g = E^^-^-^^ + E^^-^-^g. The CRHB calculations show rather 
small BCS-like pairing energies -Ebcs which slowly converge to zero with increasing spin and 
almost vanish already at / ~ 30 — 35h. Similar to Epairing, APNP(LN) significantly (by 
factor 5 — 10) increases the size of the BCS-like pairing energies. Although these energies 
decrease with increasing spin reflecting the quenching of the pairing correlations, they do 
not vanish even at highest calculated spins. 

The effective alignments ie/f between the lowest SD configurations in ^^^Hg and ^^^Pb 
calculated in CRMF, CRHB and CRHB+LN are shown in Fig. ^ and are compared with 
experiment. It is clearly seen that the CRMF and especially the CRHB results deviate 
considerably from experiment in absolute value. In addition, the considerable change of 
the slope of ie/f obtained in these calculations at Qx ~ 0.3 — 0.35 MeV, which is due to 
proton band crossing, contradicts to experimental data. Similar to the moments of inertia, 
APNP(LN) considerably improves the agreement between calculations and experiment also 
for the effective alignments. Although the values ie/f calculated in CRHB+LN deviate by 
0.4 — 0.5h from experiment, this deviation should not be considered as crucial since the 
compared bands differ by 2 protons. Moreover, the slope of ieff as a function of is rather 
well reproduced in the calculations. 



IV. THE DEPENDENCE OF THE RESULTS ON THE PARAMETRIZATION OF 
THE MEAN FIELD AND PAIRING FORCE. 

A. The dependence of the results on the parametrization of the RMF Lagrangian. 

The pairing correlations depend not only on the properties of the effective pairing force, 
but also on the single-particle level density. A full relativistic Hartree-Bogoliubov calculation 
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is therefore only meaningful if the Hartree field yields a reasonable single-particle spectrum 
28| . In order to investigate the dependence of the results on the parametrization of the 



RMF Lagrangian, the CRHB+LN calculations have been performed with the NL3 force ||62 



for the lowest SD configurations in ^^^Pb and ^^^Hg, see Figs. || and ^ In addition, the 
NLSH force has been employed in ^^^Hg, but due to slow convergence it was used only 
in a short frequency range. In all these calculations, the DIS set has been used for the 
Gogny force. 

At low rotational frequencies, total, neutron and proton kinematic and dynamic moments 
of inertia obtained in the calculations with NL3 are smaller than the corresponding quantities 
calculated with NLl, see Fig. ^. However, the increase of J*^^^ and J*^^^ as a function of 
rotational frequency is larger in the calculations with NL3 compared with the ones employing 
NLl. Thus at some frequencies, they reach each other. At even higher frequencies, the J*^^-* 
and J'-^^ values calculated with NL3 become larger than the ones calculated with NLl. It 
is also clear that the NLl force provides better agreement with experimental kinematic and 
dynamic moments of inertia than NL3. The results of the calculations with NLSH (see Fig. 
^) are in even larger disagreement with experiment than the ones with NL3. One should 
note however that the NL3 force provides better reproduction of effective alignment in the 
i94jjgy'i94p|^ pair compared with the NLl force, see Fig. 

In the NLl parametrization, the Qt values are nearly constant as a function of both 
in ^^"^Hg and ^^^Pb nuclei. The Qti^^'^Ph) is approximately 2eb larger than (5t(^^^Hg). The 
situation is different when NL3 force is used in the calculations. At f^a, ~ MeV, Qt(^^^Pb) 
is approximately equal to Qt(^^^Hg). This fact is most likely related to the = 112 SD shell 
gap which is more pronounced in the NL3 parametrization (see Fig. ^ and discussion in Sect. 
IH). However, with increasing rotational frequency Qt(^^^Pb) increases considerably with a 
maximum gain of ~ 1.5 eb at Qx ~ 0.35 MeV (Fig. ^c). On the contrary, with exception 
of the band crossing region the evolution of (5t(^^^Hg) as a function of Qx is similar to the 
one seen in the calculations with NLl. Comparing different parametrizations of the RMF 
Lagrangian, it is clear that NLl produces the largest values of Qt, while NLSH the smallest 
(Fig. ^)). This tendency has already been seen in the A ~ 60 and A ~ 140 — 150 regions 
of super deformation [|l3,|l51 • 

In the calculations with NL3, proton and neutron pairing energies are similar at low 
rotational frequencies (see Figs. ||a,b). At fix > 0.15 MeV, proton pairing energies are 
larger than neutron ones. On the other hand, in the calculations with NLl neutron pairing 
energies are larger than proton ones at all rotational frequencies. The comparison of single- 
particle energies at Qx = 0.0 MeV obtained with the NLl and NL3 parametrizations of the 
RMF theory (see Fig. 0) suggests that this is due to the larger A^ = 112 and the smaller 
Z = 80 SD shell gaps in the calculations with NL3. This leads to an additional quenching of 
neutron pairing correlations and to an increase of proton pairing correlations as compared 
with the case of the NLl force. 



B. The dependence of the results on the parametrization of the Gogny force. 

In the present section, we will study how the results of the calculations depend on the 
parametrization and the strength of the Gogny force. In all calculations given in this section, 
the set NLl is used for the RMF Lagrangian and approximate particle number projection 
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is performed by means of the Lipkin-Nogami method. Such a study is motivated by the 
fact that a precise quantitative information on the pairing correlations is not easy to extract 
in nuclei. There are no simple physical processes allowing one to isolate completely the 
pairing effects from the rest and to use them for the fit of the interaction in the particle- 
particle channel. Different sets of the Gogny force have been obtained by the fit to the 
properties of finite nuclei. As a result, these fits are more sensitive to the properties in 
the particle-hole channel than in the particle-particle channel since the pairing energies 
represent only a small portion of the total binding energies. In addition, apriori it is not 
clear that existing parametrizations of the Gogny force should be reliable in conjuction with 
the RMF theory. One should note that the moments of inertia of rotating nuclei are very 
sensitive to the properties of the pairing interaction ||2^J6^ . Considering that the rotation- 
vibration coupling is small in strongly deformed nuclei |2^;^^, one can try to use this fact 
for the definition of the best parametrization of the Gogny force which has to be used in the 
particle-particle channel in conjuction with the RMF theory. 

The parameter set Dl has been defined in Ref. |^5|] based on the features of the ^^O, 
^^Ca, ^°Zr nuclei and the pairing properties of the Sn isotopes. It turns out that in non- 
relativistic calculations with the Gogny force this set overestimates the strength of the pairing 
correlations. The Dl' set ||6^ differs from the Dl set only in the strength of the spin-orbit 
interaction and thus it will not be employed in the present calculations since only the central 
force with two-finite range Gaussian terms (see Eq. (|l^)) of the original Gogny force 
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used in the CRHB theory. The DIS set |34|J67[] differs, as follows from non-relativistic studies, 
from Dl by improved surface properties and by producing smaller pairing correlations. It 
is used in most of the calculations in the present manuscript. Recently another set of the 
Gogny force has been suggested in Ref. So far it has not been applied for detailed 



studies of finite nuclei in non-relativistic approaches. However, comparing it with Dl and 
DIS (see Table |), one can conclude that it has larger similarities with Dl than with DIS. 

The dynamic and kinematic moments of inertia of the lowest SD configuration in ^^^Hg 
calculated with Dl, DIS and DIP sets are compared with experiment in Fig. |1^. One can 
see that only the set DIS provides good agreement with experimental data. The values of 
J^^^ calculated with Dl and DIP are appreciable below both the experiment and the results 
obtained with DIS, see Fig. |lO|b. This is mainly due to smaller J^^^ values for neutrons. 
At low rotational frequencies, the J^^^ values obtained with Dl and DIP are lower than 
both experimental data and the ones obtained with DIS. However, the increase of J'-^-' as a 
function of rotational frequency is larger in the calculations with Dl and DIP and thus at 
fix ~ 0.32 MeV they become larger than both experiment and the values calculated with 
DIS. In the neutron band crossing region at ^l^ ~ 0.42 MeV, the J*^^-* values calculated with 
Dl and DIP are larger than both the ones obtained with DIS and the experiment. It is 
interesting to note that all three sets give very similar neutron band crossing frequencies. 
Similar to the case of J^^\ the differences between J^'^^'s calculated with Dl and DIP on 
the one side and with DIS on the other side are connected mainly with different alignment 
patterns in neutron subsystem. 

Calculated neutron and proton pairing energies ^pai^j^g are displayed in Fig. |ll|b. The re- 
sults of the calculations with different sets show a similar behavior as a function of rotational 
frequency, but the absolute values of pairing energies strongly depend on the parametriza- 
tion of the Gogny force. The strongest pairing correlations are provided by the Dl set, while 
the smallest ones are obtained with DIS. The pairing energies calculated with DIP are in 
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between the ones obtained with DIS and Dl. Comparing transition quadrupole moments 
Qt obtained with different sets (see Fig. [TT|a), one can see that (i) the Dl and DIP sets 
provide very similar values of Qt and (ii) at high rotational frequencies the values of Qt 
calculated with all three sets come very close to each other reflecting the decrease of pairing 
correlations. 

An alternative way to study the role of pairing correlations for different physical observ- 
ables is to look how the results of the calculations are affected when the strength of the 
Gogny force is changed. This is done by introducing the scaling factor / into the Gogny 
force (see Eq. (plTf)). Such an investigation is in part motivated by the fact that in many 
studies performed within the relativistic Hartree-Bogoliubov theory with Gogny forces in 
the pairing channel and with no particle number projection, the strength of the Gogny force 



is increased by the factor 1.15, see for example Ref. The results of the CRHB calcula- 
tions for the lowest SD configuration in ^^^Hg with the scaling factors / = 0.9, / = 1.0 and 
/ = 1.1 for the DIS Gogny force are shown in Figs. |12] and [T^. The effect of the scaling of 
the Gogny force is especially drastic on the pairing and rotational properties of the configu- 
ration under investigation (see Figs. and |r2|). The increase (decrease) of the strength of 
the Gogny force by 10% leads to approximately twofold increase (decrease) of the absolute 
values of proton and neutron pairing energies, see Fig. |13]b. On the contrary, the increase 
(decrease) of the strength of the Gogny force leads to a significant decrease (increase) of the 
kinematic moment of inertia (Fig. 0b). One should note that the impact of the change of 
the strength of the Gogny force on the proton and neutron J^^^ values is different. The effect 
of the scaling of the Gogny force on the dynamic moment of inertia is more complicated. 
While at low rotational frequencies it is similar to the one seen for the kinematic moment of 
inertia, it is completely different in the paired band crossing regions as seen in the proton, 
neutron and the total dynamic moments of inertia (Fig. |l2|a). Contrary to the low frequency 
range, stronger pairing correlations lead to larger J^^^ values. In addition, as it is seen on 
the example of the neutron subsystem, stronger pairing correlations (a larger strength of the 
Gogny force) lead to the shift of the neutron paired band crossing to higher frequencies and 
make it sharper. A larger strength of the Gogny force leads also to smaller values of Qt at 
low fia; and to larger Qt at high fi^;, see Fig. p^a. 

In addition, different strengths and different sets of the Gogny force lead to different 
quasiparticle spectra in the vicinity of the Fermi level. In part, this effect is caused by the 
different calculated equilibrium deformations. It is reasonable to expect that this dependence 
of the quasiparticle spectra upon the parametrization and the strength of the Gogny force 
will have a pronounced impact on the band crossing frequencies in the SD configurations in 
odd and odd-odd nuclei. 



C. Concluding remarks. 

The difference between the transition quadrupole moments Qt calculated with different 
sets of the RMF forces using the same set of the Gogny force is appreciable (Figs. ||c,d). On 
the contrary, the difference between QtS obtained with three different sets (Fig. |ll|) or with 
three different scalings of the strength (Fig. |T3|a) of the Gogny force is smaller reflecting 
the fact that the equilibrium deformations are mostly defined by the properties of effec- 
tive interaction in the particle-hole channel. The experimental uncertainties on transition 
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quadrupole moments (see discussion in Ref. [^) prevent the use of calculated Qts for the 



selection of a better parametrization of the RMF Lagrangian for the particle-hole channel 
and the Gogny force for the particle-particle channel. Thus this selection can be based only 
on the rotational properties of the SD bands under study such as the kinematic and dynamic 
moments of inertia which are very accurately defined in experiment. 

Comparing different parametrizations of the Gogny force, it is clear that the Dl and DIP 
sets in connection with the NLl force provide too strong pairing correlations and give very 
similar results for kinematic and dynamic moments of inertia, which deviate appreciable from 
experiment. The results of the calculations with the NL3 and NLSH forces in conjuction 
with the DIS set of the Gogny force also lead to considerable deviations from experiment 
for and J^'^\ It is then reasonable to expect that the use of the Dl and DIP sets of 
the Gogny force in conjuction with NL3 and NLSH will lead to even larger deviations from 
experiment for kinematic and dynamic moments of inertia. Thus one can conclude (see also 
Ref. ||2^) that only the NLl force in conjuction with the DIS set of the Gogny force and 



APNP(LN) will lead to a reasonable description of rotational properties of SD bands in the 
A ~ 190 mass region. 



V. PROPERTIES OF YRAST SD BANDS IN EVEN-EVEN NUCLEI 

It is our believe that the strong and the weak points of a theoretical approach can be 
determined only by a systematic comparison between experiment and theory. In the present 
Section the results of a systematic investigation of yrast SD bands in even-even nuclei of the 
A ~ 190 mass region will be presented. This investigation covers all even-even nuclei in this 
region in which SD bands have been observed so far, namely, I90'i92,i94jjg^ I92,i94,i96,i98p|j 

and ^^'^Po. First, the results for Hg and Pb isotopes with N = 110, 112 and 114 will be 
presented in detail. Partial results for these nuclei have been already presented in Ref. ||20|| . 
Then the results for ^^'^Pb and ^^^Po will be discussed for each nucleus separately. All the 
calculations presented in this Section have been performed with the force NLl for the RMF 
Lagrangian and the set DIS for the Gogny force in the pp-channel and approximate particle 
number projection by means of the Lipkin-Nogami method (APNP(LN)) has been used. 



A. The TV = 110, 112 and 114 Hg and Pb nuclei 

The calculated total, neutron and proton dynamic and kinematic moments of inertia 
are shown in Figs. |1^, |15| and compared with experiment. One can see that a very suc- 
cessful description of dynamic moments of inertia of experimental bands is obtained in the 
calculations without adjustable parameters (Fig. [1^). When comparing calculated and ex- 
perimental kinematic moments of inertia, one should keep in mind that only yrast SD bands 



in ^^^Pb and ^^^Hg are definitely linked to the low-spin level scheme [l70| , |7ll , [73| . In addition, 
there is a tentative linking of the SD band in ^^^Pb |]73|. On the contrary, at present the 
yrast SD bands in i^O'i^^gg g^j^^ ^^^Pb are not linked to the low-spin level scheme yet. Thus 
some spin values consistent with the signature of the calculated lowest SD configuration 
should be assumed for the experimental bands when a comparison is made with respect 
of the kinematic moment of inertia J^^\ Taking into account that kinematic moments of 
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inertia of linked SD bands in ^^^Pb and ^^^Hg and tentatively linked SD band in ^^^Pb are 
well described in the calculations, the spin values of unlinked bands can be obtained by com- 
paring the calculated values of J*-^-* with experimental ones under different spin assignments. 
Such a comparison (see Ref. [^] for detailed figures) leads to the spin values Iq for the 
lowest states in SD bands listed in Table fT|. Under these spin assignments, the CRHB+LN 
calculations rather well describe 'experimental' kinematic moments of inertia of SD bands 
in 192,196 p]^ Hg (Fig. 0). Alternative spin assignments, which are different from the 

ones given in Table |T| by ±2h, can be ruled out since they (open circles in Fig. |l^ lead to 
considerable deviations from the results of the calculations. 

The increase of kinematic and dynamic moments of inertia in this mass region can be 
understood in the framework of the CRHB+LN theory as emerging predominantly from a 
combination of three effects: the gradual alignment of a pair of ji5/2 neutrons, the alignment 
of a pair of zi3/2 protons at a somewhat higher frequency, and decreasing pairing correlations 
with increasing rotational frequency, see also Sect. |T|. The interplay of alignments of neutron 
and proton pairs is more clearly seen in Pb isotopes where the calculated J*^^-* values show 
either a small peak (for example, at ~ 0.44 MeV in ^^^Pb, see Fig. or a plateau 
(at ilx ~ 0.4 MeV in ^^^Pb, see Fig. With increasing rotational frequency, the J'-^-' 
values determined by the alignment in the neutron subsystem decrease. The maximum in 
J(^) of the neutron subsystem is reached at fi^ ^ 0.44, 0.41, 0.39, 0.425 and 0.39 MeV in 
I92,i94,i96p|^ and i94,i96jjg^ respectively. The decrease in these frequencies with increasing A'" 
within each isotope chain correlates with the decrease of the transition quadrupole moments 
Qt (Fig. 4 in Ref. [^). However, the fact that the maximum in neutron J^'^^ is obtained at 
the same frequencies in Hg and Pb isotones, which have Qt values differing by 1 — 1.5 eb 
(Fig. 4 in Ref. indicates that the alignment process is very complicated and depends 

not only on the equilibrium deformation but also on the position of the Fermi level. 

The decrease of neutron J^^^ at frequencies higher than ^l^ ~ 0.4 MeV is in part compen- 
sated by the increase of proton J*^^-* due to the alignment of the ^13/2 proton pair. This leads 
to the increase of the total J'-^-'-value at Vlx > 0.45 MeV in the Pb isotopes, while no such 
increase has been found in the calculations after the peak up to Qx = 0.5 MeV in ^^^Hg, see 
Fig. |l^. Thus one can conclude that the shape of the peak (plateau) in total J*^^-* in the band 
crossing region is determined by a delicate balance between alignments in the proton and 
neutron subsystems which depends on deformation, rotational frequency and Fermi energy. 
It is also of interest to mention that the sharp increase in J(2) of the yrast SD band in ^^^Hg 
is also reproduced in the present calculations. In the calculations, this increase is due to a 
two-quasiparticle alignment associated with the z/[761]3/2 orbital. One should note that the 
calculations slightly overestimate the magnitude of J'-^-* at the highest observed frequencies. 
Possible reasons could be the deficiencies either of the Lipkin-Nogami method fSlI] or of the 
cranking model in the band crossing region or both of them. In addition, the calculations 
do not reproduce the sudden decrease in J^"^^ at the bottom of the SD band in ^^^Pb, the 



origin of which is not understood so far, see Ref. [|75 |. 

In order to compare relative properties of the calculated and experimental dynamic mo- 
ments of inertia in more detail we tilt them by extracting the frequency dependent term 
from J^^-*. The resulting quantities (J*-^-* — ISOilx) are shown in Fig. |16| for Hg (top panel) 
and Pb (middle panel) isotopes as well as for = 112 (bottom panel) isotones. Since the 
difference of the dynamic moments of inertia of two bands is proportional to the derivative 
of the effective (relative) alignment ieff of these bands (see Refs. [p2| , p!8| ) for details), the ef- 
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fective alignments of compared bands will also be presented here. The relative properties of 
the dynamic moments of inertia of the yrast SD bands in i92,i94jjg nuclei are very well repro- 
duced in the calculations. Both in calculations and in experiment, J*^^)(^^^Hg):^J*^^)(^^^Hg) 
at Vt^ < 0.3 MeV, while at higher frequencies j(2)(i94Hg)> j{2)(i92Hg) (Figs. |l6|a,b). This 
result correlates with the fact that effective alignment ig// in the ^^^Hg/^^^Hg pair and es- 
pecially its slope is also reproduced rather well in the calculations (see Fig. p!7|d). On the 
contrary, while the properties of J*^^)(^^°Hg) relative to J*^^''(^^^Hg) are reasonably well re- 
produced at frequencies fl^ > 0.3 MeV, the situation is different at lower frequencies. There 
the difference between J^^-* values stays almost constant around 6 MeV~^ in experiment, 
while it is decreasing to zero at Qx ~ 0.2 MeV in the calculations (see Figs. |TB|a,b). This 
feature correlates with the fact that the effective alignment in this pair of SD bands is not 
well reproduced in the calculations (see Fig. 0b). 

The dynamic moments of inertia of yrast SD bands in i92,i94p]^ almost identical 
in experiment. This feature (Figs. |l6|c,d) and the effective alignment in the ^^^Pb/^^^Pb 
pair (Figs. |T7|a) are rather well reproduced in the calculations. In experiment, the dynamic 
moments of inertia of SD bands in i94,i96p]^ ^^^^ almost identical at ~ 0.1 MeV. With 
mcreasmg 

Qx, J(^n^^^Pb) decreases below J^^^^^^Ph) with a maximum difference between 
them of ~ 5 MeV~^ being reached at fl^ ~ 0.22 MeV and then this difference becomes 
smaller up to ~ 0.32 MeV where the dynamic moments of inertia of both bands coincide 
(Fig. ^c). At even higher frequencies, J'^^)(^^^Pb)> J*^^)(^^^Pb). The calculations reasonably 
well reproduce the general features, however, somewhat underestimate the difference between 
J*^^)'s of these two bands at medium frequencies and overestimate this difference at highest 
observed frequencies. The experimental effective alignment in the ^^^Pb/^^^Pb pair is very 
well reproduced at low frequencies, while the difference between experiment and calculations 
is somewhat larger at high frequencies (Fig. [T7|c). 

The experimental dynamic moments of inertia of SD bands in ^^^Pb and ^^^Hg are iden- 
tical at < 0.2 MeV, while at higher frequencies the condition j(^)(^^^Pb)>J*^^)(^^^Hg) 
holds (Fig. |16|e). The relative properties of dynamic moments of inertia of these two bands 
at > 0.2 MeV are rather well reproduced in the calculations (Fig. 0f). At lower frequen- 
cies, there is however some difference between the calculated J*^^-* values for these two bands 
of ~ 2 MeV~^. The effective alignment in the ^^^Hg/^^^Pb pair and especially its slope as a 
function of fix is well reproduced in the calculations (Fig. [l7|e). 

The effective alignments between other pairs of SD bands which have not been discussed 
before are also shown in Fig. p!7|g,f,i. It is seen that the effective alignment in the ^^^Hg/^^^Pb 
pair is reasonably well reproduced in the calculations, although they somewhat overestimate 
the absolute value of ieff- Since the spins of these two bands are determined experimentally, 
this result can be considered as a first direct justification of the reliability of the effective 
alignment approach used frequently for the configuration and spin assignments of SD bands 
in different mass regions, see Refs. JTH|,|T2[ for details. Comparing the experimental and 
effective alignments (Fig. |17[) and taking into account that compared bands differ by at 
least 2 particles, one can conclude that considerable deviations from experiment are seen 
only in the cases of the effective alignments in the ^^°Hg/^^^Pb, ^^^Hg/^^^Pb and ^^°Hg/^^^Hg 
pairs (Figs. |l^g,f,b). In order to find are observed deviations from experiment related to the 
particle-hole or the particle-particle channel of the CRHB theory the investigation of SD 
bands in odd nuclei of this mass region are needed. Such an investigation is in progress and 
its results will be reported later. The investigation of relative alignments of SD bands in 
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this mass region has been performed in non-relativistic approaches such as the total routhian 
surface Strutinsky-type approach based on the Woods-Saxon potential and Skyrme-Hartree- 



Fock-Bogoliubov approach in Ref. ||83[. Similar to our case, some discrepancies between 
theory and experiment have been found. 

The calculated neutron and proton pairing energies E^^^^^^g are summarized in Fig. |T^. 
In all nuclei, they decrease with increasing rotational frequency reflecting the quenching of 
pairing correlations due to the Coriolis antipairing effect. In addition, neutron pairing is 
stronger than proton pairing. Neutron pairing energies in i^O'is^pjg almost coincide up to the 
band crossing seen in ^^°Hg. A similar situation exists also in i92,i94p|^ where the difference 
between calculated neutron pairing energies does not exceed 0.2 MeV. These energies are 
smaller than the ones in i^O'^s^j^g ]^-y ^ g^^ MeV at VL^ = 0.0 MeV, while at high rotational 
frequencies neutron pairing energies are similar in these Pb and Hg nuclei. On the contrary, 
neutron pairing energies in the = 114 Hg and Pb nuclei are larger than the ones in the 
nuclei with = 112 by 0.5-0.8 MeV dependent on the rotational frequency and they do 
coincide at fi^. > 0.3 MeV. A similar trend is seen also for proton pairing energies where, 
however, the difference between the pairing energies in different nuclei is smaller than in the 
case of neutron pairing energies and it does not exceed 0.5 MeV. 



B. The nucleus ^^^Pb 
The results of the calculations for the lowest SD band in ^^'^Pb are shown in Figs. ^ 



and The calculated kinematic moment of inertia agrees reasonably well with experiment 
up to ~ 0.32 MeV, while considerable disagreement is seen at higher frequencies. The 
increase of calculated J^^^ as a function of fl^ is larger in ^^^Pb than in ^^^Pb (Fig. p!9|b). 
As a result, the calculated dynamic moment of inertia in ^^'^Pb is larger than the one in 
^^^Pb at Qx > 0.2 MeV and smaller at < 0.2 MeV. On the contrary, the experimental 
data show the opposite trend for the dynamic moment of inertia of the ^^®Pb band having a 
much smaller increase as a function of rotational frequency (Fig. p^a). Thus the results of 
the calculations for ^^^Pb do not reproduce neither absolute rotational properties of SD band 
nor their relative properties with respect to SD band in ^^^Pb. A similar problem with the 
reproduction of the properties of the ^^^Pb band exists also in the cranked Nilsson-Strutinsky 



Lipkin-Nogami calculations presented in Refs. p6| , |79| . The investigation of neighbouring odd 



nuclei within the CRHB+LN theory is needed in order to understand better the origin of 
these problems. The calculated values of Qt for SD band in ^^^Pb are by ~ 1 eb smaller 
than in ^^^Pb (see Fig. pO|a) which is in agreement with the decrease of Qt with increasing A^ 
seen in the lighter Pb isotopes (see Fig. 4 in Ref. |^). Precise measurements of the relative 
transition quadrupole moments of the yrast SD bands in i96,i98p|^ nuclei and a theoretical 
study of neighboring odd nuclei can be useful for the understanding of the present problems 
with the description of J^^^ and J^'^\ Pairing energies in ^^^Pb are somewhat larger than in 
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Pb (Fig. 20 b) which agrees with the trend of the increase of pairing correlations within 



the isotopic chain with increasing neutron number A^ (see Section |V A|) . 
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C. The nucleus ^^^Po 

The CRHB+LN calculations for the lowest SD configuration in ^^^Po very well describe 
experimental dynamic and kinematic moments of inertia of the yrast SD band in this nucleus 
(Fig. pT|a). At Qx ~ 0.34 MeV, the calculations predict sharp increase in dynamic moment of 
inertia caused mainly by the alignment of the lowest proton hyperintruder vr [770] 1/2 orbital 



(see Fig. |22|). The absolute values of proton and neutron pairing energies and their behavior 
as a function of rotational frequency (Fig. pl|b) are similar to the ones seen in other nuclei. A 
specific feature of this nucleus, which has not been observed in other nuclei, is considerable 
increase of the transition quadrupole moment Qt (by ~ 3 eb) with increasing rotational 
frequency (Fig. ^). At = 0.0 MeV, the Qt values in ^^^Po are larger than the ones in 
isotonic ^^^Pb by ~ 1.0 eb (see Fig. |2T| c and Fig. 4 in Ref. |2^). One should note that the 
results of the GCM+GOA calculations based on the Gogny force also show the same feature 
|60| . Finally, the experimental effective alignment in the ^^^Pb/^^'^Po pair is reasonably well 



described in the calculations (Fig. ^ 



D. Mass hexadecupole moments Qt 

Experimental information on mass hexadecupole moments Q40 is not available so far for 
SD bands in any mass region. Thus only the results of the calculations for this quantity 
are presented in Fig. These results have been obtained with the NLl force for the RMF 
Lagrangian and the set DIS for the Gogny force if it is not specified otherwise. Note that we 
do not show the results of the calculations obtained with either different parametrizations or 
different scalings of the Gogny force. The Q40 values calculated with no pairing (dotted lines 
in Figs. [23|a,b) show a gradual decrease with an increase of rotational frequency. A similar 
trend is also seen in the results of the calculations with no APNP(LN) (solid lines with solid 
squares in Figs. |23|a,b). The sharp change of the slope of the Q^oi^x) curve seen in ^^^Pb at 
Qx ~ 0.35 MeV (Fig. p3|b) is related to a sharp band crossing associated with the collapse of 
proton pairing correlations (see Sect. |T|). The results of the calculations with APNP(LN) 
(open symbols in Fig. show a somewhat different trend. With the exception of ^^'^Po, the 



calculated Q40 values stay nearly constant or smoothly increase with increasing rotational 
frequency up to Qx ~ 0.35 MeV. Above this frequency they decrease with increasing Qx- In 
^^^Po, the Q40 values increase with increasing Qx in the whole calculated frequency range 
(Fig. ^3]c). This increase is especially pronounced in the band crossing region. Within the 
isotopic chain the increase of neutron number N leads to the decrease of Q40. The Q^q 
values increase within isotonic chain with increasing proton number Z. The results of the 
calculations with NL3 and NLSH lead to smaller values of Q40 compared with the ones 
obtained with NLl force (Fig. p3|a,b). The features discussed above are very similar to 
the ones obtained for the transition quadrupole moment Qt, see the discussion in previous 
sections. 



E. Particle number fluctuation ((AA^)^ 



The basis assumption behind the Kamlah expansion to second order used in the 
derivation of the Lipkin-Nogami method is that the system is well pair-correlated which 
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means that the particle number fluctuation in the unprojected wave function ((AiV)^) is 
large. These quantities obtained in the CRHB+LN calculations with the NLl and NL3 
forces for the RMF Lagrangian and the DIS set for the Gogny force are shown in Fig. 
for all nuclei studied in the present manuscript. The particle number fluctuations de- 



crease more or less smoothly with increasing rotational frequency indicating the quenching 
of pairing correlations due to the Coriolis antipairing effect. One should note that even at 
the highest rotational frequencies these fluctuations remain reasonably large thus indicating 
that the approximate particle number projection by means of the Lipkin-Nogami method 
still remains within the applicability range of the Kamlah expansion. The values of ((AA^)^) 
in neutron and proton subsystems correlate with the pairing energies -Epairmi; calculated in 
these subsystems. For example, in the calculations with the NLl force the particle number 
fluctuations are larger for neutrons than for protons (see Fig. ^4]) which correlates with the 



fact that the absolute values of pairing energies are larger for neutrons (see Figs. |20[ b 
and |21[b). The situation is somewhat different in the calculations with the NL3 force, where 
at medium and high rotational frequencies the proton subsystem is more pair-correlated 
than the neutron one as reflected in the particle number fluctuations (Fig. and pairing 
energies (Fig. Ha,b). 



VI. CONCLUSIONS 

The formalism of the Cranked Relativistic Hartree-Bogoliubov theory with and without 
approximate particle number projection before variation by means of the Lipkin-Nogami 
method is presented in detail. The relativistic mean field theory is used in the particle-hole 
channel of this theory, while a non- relativistic finite range two-body force of Gogny type is 
employed in the particle-particle (pairing) channel. Considering that the pairing is a genuine 
non-relativistic effect which plays a role only in the vicinity of the Fermi surface, the use of 
the best non-relativistic force in the pairing channel seems well justified. 

Its applicability to the description of rotating nuclei and the main features of this theory 
have been studied on the example of the yrast superdeformed bands observed in even-even 
nuclei of the A ~ 190 mass region. The main conclusions emerging from this study are the 
following: 

(i) The calculations without particle number projection do provide only a poor descrip- 
tion of experimental rotational features such as the kinematic J*^^-* and the dynamic J*-^-* 
moments of inertia. The calculated kinematic moments of inertia are larger than the ex- 
perimental values. The same is also true for the dynamic moments of inertia at low and 
medium rotational frequencies. The calculations without particle number projection lead 
to a unphysical collapse of pairing correlations as has been seen in the proton subsystem of 
^^^Hg and ^^^Pb. As was shown by subsequent calculations with particle number projection 
these problems are related to the poor treatment of pairing correlations. 

(ii) Approximate particle number projection by means of the Lipkin-Nogami method 
considerably improves an agreement with experiment. The correlations induced by the 
Lipkin-Nogami method produce the desired effects, namely, (a) to increase the strength of 
the pairing correlations, (b) to diminish the kinematic moments of inertia at all frequencies, 
(c) to diminish the dynamic moments of inertia at low and medium rotational frequencies 
and (d) to delay proton and neutron alignments to higher frequencies. In addition, there is 
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no collapse of pairing correlations in the whole rotational frequency range under investiga- 
tion. Systematic calculations with the NLl force for the RMF Lagrangian and the DIS set 
for the Gogny force have been performed for all even-even nuclei in which SD bands have 
been observed so far. With an exception of ^^^Pb, an excellent description of the rotational 
properties of yrast SD bands such as the dynamic and kinematic moments of inertia is ob- 
tained in a way free from adjustable parameters. It was concluded that the investigation of 
the structure of SD bands in neighboring odd nuclei is needed in order to better understand 
the problems seen in this nucleus. 

(iii) It has been investigated how much the results of calculations with approximate par- 
ticle number projection by means of the Lipkin-Nogami method depend on the parametriza- 
tion of the RMF Lagrangian and the Gogny force. It was found that the combination of the 
NLl set for the RMF Lagrangian and the DIS set for the Gogny force produce very good 
agreement with experimental rotational properties. The Dl and DIP sets of the parameters 
for the Gogny force produce too strong pairing correlations and thus fail to describe properly 
the rotational properties of SD bands. An unexpected result of the present investigation is 
the fact that the NL3 force, which is believed to be the best RMF force for the description of 
nuclear properties far from beta-stability region, provides a less accurate description of the 
rotational properties of SD bands in the A ~ 190 mass region compared with the force NLl. 
This is most hkely related to the single-particle spectra produced by this force in the SD 
minimum. It is possible that the improvement of the description of the isospin properties 
far from beta-stability region obtained in the NL3 force as compared with the NLl force is 
reached at the cost of some worsening of the description of the single-particle spectra close 
to beta-stability region. 

(iv) The dependence of the results on the strength of the Gogny force in the pairing 
channel has been studied in the calculations with approximate particle number projection 
by means of the Lipkin-Nogami method. It was found that a change of the strength by 
±10% has significant impact on the rotational properties of SD bands such as the kinematic 
and dynamic moment of inertia. Contrary to previous investigations within the Relativistic 
Hartree-Bogoliubov theory at no rotation and with no particle number projection, where the 
strength of the Gogny force has been increased by factor 1.15, here a very good description 
of rotational properties has been obtained with no modification of the strength of the Gogny 
force. 

(v) The results of the calculations with particle number projection indicate the general 
trend of a decrease of the average transition quadrupole {Qt) and mass hexadecupole {Q40) 
moments in the isotopic chain with increasing neutron number N and an increase of these 
quantities within the isotonic chain with increasing proton number Z. 

(vi) Neutron and proton pairing energies in all the calculated nuclei decrease with in- 
creasing rotational frequency reflecting the quenching of pairing correlations due to the 
Coriolis anti-pairing effect. The pairing energies Epamng calculated in the NLI+DIS+LN 
scheme do not show considerable variations as a function of Z and N\ the maximum differ- 
ence between pairing energies calculated in two different nuclei is around 16%. A smooth 
increase of absolute values of pairing energies with increasing N is observed in the isotopic 
Pb and Hg chains which correlates with the decrease of the transition quadrupole moments 
Qt- The size of the pairing energies strongly depends on the parametrization of the Gogny 
force and the RMF Lagrangian, as well as on the strength of the Gogny force. 

(vii) The difference between the experimental dynamic moments of inertia of two bands 
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and the effective alignment ig// between these bands is reasonably well reproduced in most 
of the cases. Further investigation of neighboring odd nuclei is needed in order to find the 
origin of the remaining discrepancies. 

Finally, the present work should be considered as one of the first steps in the investigation 
of rotating nuclei in the pairing regime within the framework of the RMF theory. Different 
tasks definitely lie ahead. For example, an investigation of the rotational bands based on 
one- and multi-quasiparticle configurations in the A ~ 190 mass region is mandatory in 
order to see how the present theory can reproduce the effects connected with the blocking 
of one or several quasiparticle orbitals. One can expect that the Lipkin-Nogami method is 
a reasonably good approximation to the exact particle number projection in the regimes of 
strong pairing correlations as it holds in the case of the A ~ 190 mass region. It remains 
to investigate if this method is also a good approximation in the regimes of weak pairing 
correlations typical at high rotational frequencies in SD bands of other regions of the periodic 
table, such as the A ~ 60, A ~ 130 and A ~ 150 mass regions. 
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FIG. 1. Neutron (top) and proton (bottom) single-particle energies (routhians) in the 
self-consistent rotating potential as a function of rotational frequency Cl^. They are given along 
the deformation path of the lowest SD configurations in ^^^Hg (left panels) and ^^^Pb (right pan- 
els) obtained in the calculations neglecting pairing and with the NLl parametrization of the RMF 
Lagrangian. The calculated transition quadrupole moments Qt are shown on the top of the figure 
as {Qt = Qti^x = 0.0) Qti^x = 0.5)). The notation of lines used for routhians is also indi- 
cated on the top of the figure. Chemical potentials A (Aln) obtained in the paired calculations 
without (with) particle number projection by means of the Lipkin-Nogami method are shown by 
long-dashed lines. 
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mations of the lowest SD configuration in ^^^Hg. The single-particle orbitals are labeled by means 

of the asymptotic quantum numbers [A^n^A]$7 (Nilsson quantum numbers) of the dominant com- 
ponent of the wave function. Solid and dashed lines are used for the positive and negative parity 
states, respectively. The calculated charge quadrupole Qq and mass hcxadccupolc Q^^ moments 
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on the pairing and particle number projection illustrated on the example of the lowest SD con- 
figurations in ^^^Pb and ^^^Hg. Experimental data are shown by solid unlinked symbols. The 
calculations have been performed with no pairing (indicated as NLl [CRMF]), without particle 

number projection (NLI+DIS) and with approximate particle number projection by means of the 
Lipkin-Nogami method (NLI+DIS+LN). Proton and neutron contributions to the kinematic and 
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FIG. 7. The dependence of the effective ahgnment ig// on the pairing and the particle number 
projection (panel (a)) and on the parametrization of the RMF Lagrangian (panel (b)). In both 
panels, the experimental data are shown by unlinked solid circles, while the results of the calcu- 
lations by lines of different type. The experimental effective alignment between bands A and B is 
indicated as "A/B". Band A in the lighter nucleus is taken as a reference, so the effective alignment 
measures the effect of the additional particle(s). In panel (a), the results of the calculations with no 
pairing (indicated as NLl [CRMF]), without particle number projection (NLI+DIS) and with ap- 
proximate particle number projection by means of the Lipkin-Nogami method (NLI+DIS+LN) are 
compared with the experimental effective alignment in the ^^^Hg/^^^Pb pair. Panel (b) compares 
effective alignments in the ^^^Hg/^^^Pb pair obtained with NLl and NL3 forces with experimental 
data. 
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FIG. 8. The dependence of the kinematic (J*-^^) and dynamic (J^^^) moments of inertia of the 
lowest SD configurations in ^^^Pb and ^^^Hg on the parametrization of the RMF Lagrangian. The 
calculations have been performed with the sets NLl, NL3 and NLSH. In all calculations, the DIS set 
of parameters has been used for the Gogny force and approximate particle number projection has 
been performed by means of the Lipkin-Nogami method. Due to slow convergence the calculations 
with the NLSH set have been done only in the case of ^^^Hg and only in the rotational frequency 
range 0.15 — 0.25 MeV. In this frequency range the dynamic moment of inertia calculated with 
NLSH merge together with the ones obtained with NL3 and NLl and thus it is not shown. No 
convergence is obtained in the calculations with NL3 above > 0.43 MeV in the case of ^^^Pb. 
Note that in the case of ^^^Hg, the calculations with the NL3 force considerably overestimate the 
experimental dynamic moment of inertia in the band crossing region: the calculated value of J^^^ 
reaches 202 MeV~^ at 0.445 MeV. In panel (c), the proton contribution to J^^^ is shown by a 
short-dashed line in order to make it more visible. 
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FIG. 10. The dependence of the kinematic and dynamic moments of inertia of the lowest SD 
configuration in -^^^Hg on the parametrization of the Gogny force. The results with the sets Dl, 
DIP and DIS for the Gogny force are shown. 
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FIG. 11. The dependence of the transition quadrupole moment Qt and neutron and proton 
pairing energies Epairing of the lowest SD configuration in ^^^Hg on the parametrization of the 
Gogny force. The results with the sets Dl, DIP and DIS for the Gogny force are shown. 
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FIG. 12. Calculated and experimental dynamic and kinematic moments of inertia for yrast SD 
band in ^^^Hg. Experimental data are shown by solid unlinked symbols. The calculations have 
been performed with scaling of the strength of the Gogny force: scaling factors 0.9, 1.0 and 1.1 
have been used. Neutron and proton contributions into kinematic and dynamic moments of inertia 
are indicated by the letters v and tt, respectively 
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FIG. 13. Transition quadrupole moments Qt and neutron and proton pairing energies Epairing 
of the lowest SD configuration in ^^^Hg calculated with different scaling of the strength of the 
Gogny force. The notation of lines is given in the figure. 
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FIG. 14. Dynamic moments of inertia J^^^ of observed (solid circles) yrast SD bands in the 
190,192,194Hg and 192,194,196pb nuclei versus the ones of calculated lowest in energy SD configura- 
tions. Solid lines show the total calculated dynamic moments of inertia J^^-*, while long-dashed 
and dash-dotted lines show the contribution in J^^^ from neutron and proton subsystems. The 
experimental data are taken from Refs. @ (^^^Pb), C^^Fh), [|| {^^^Ph), ^ (^^°Hg), 

(i92Hg) and (is^Hg). 
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FIG. 15. The same as Fig. 14 but for kinematic moments of inertia J^^-*. The experimental 
moments of inertia of hnked or tentatively linked SD bands are shown by solid circles. In other 
cases (^^^Pb, i^O-is^Hg), the 'experimental' kinematic moments of inertia J*-^^ are shown for three 
different spin values of the lowest state Iq in the SD band, see text for details. The values being in 
best agreement with the calculations are indicated by solid circles, while open circles are used for 
alternatives. 
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FIG. 16. Comparison of experimental and calculated dynamic moments of inertia J^"^^ in the 
Hg (top panels) and the Pb (middle panels) isotopes as well as in the N = 112 isotones (bottom 
panels). The same type of symbols is used for experimental bands (left panels) and their theoretical 
counterparts (right panels) . In order to show the differences between different bands in more detail, 
the frequency dependent term is extracted from J^^^ , see text for details. 
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FIG. 18. Calculated neutron and proton pairing energies for the lowest SD configurations in 
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FIG. 19. Calculated and experimental dynamic and kinematic moments of inertia for yrast SD 
band in ^^^Pb. The results of the calculations are shown by solid lines. Experimental data for 
iQSpij )-,a,nd is taken from Ref. |79| and 'experimental' kinematic moments of inertia are shown for 4 
different spin values for lowest state Iq. The J*-^^ values being in best agreement with results of the 
calculations at low spin are shown by solid circles. For comparison, the same quantities obtained 
for the lowest SD configuration in ^^^Pb are shown by dashed lines. The experimental values of 
j(2) for yrast SD band in ^^^Pb are shown by open squares. 



53 



21 



O 19 - 



18 



' 1 


' 1 ' 1 ' 1 ' _ 
(a) - 

s 






1 1 


.1,1.1,^" 



> 







-5 



- F -10 



TO 
Q. 



LU 



-15 



. ' 1 ' 1 


' 1 ' 


1 1 


- '^'^^ 




xy 
y / 










_ ^ ^ 

■ , 1 , 1 


, 1 . 


1 1 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

Rotational frequency [MeV] 

FIG. 20. Calculated transition quadrupolc moments Qt and proton and neutron pairing energies 



of the lowest SD configuration in ^'^^Pb. For comparison, the same quantities obtained for 
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the lowest SD configuration in ^^^Pb are shown by dashed lines. 
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FIG. 21. The results of the calculations for the lowest SD configuration in ^^^Po. Dynamic and 
kinematic moments of inertia, proton and neutron pairing energies Epairing, transition quadrupole 
moment Qt and effective alignment ie// in the ^^^Fh/^^^Po pair are shown in panels (a), (b), (c) 
and (d) respectively. The experimental data for the kinematic and dynamic moments of inertia 
are taken from Ref. and shown by unlinked open and solid circles in panel (a). In panel (d), 
the experimental effective alignment is shown by unlinked solid circles. 
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FIG. 24. Particle number fluctuations ((iV)^) as a function of rotational frequency. Solid and 
dashed lines are used for neutrons and protons, respectively. Thick and thin lines are used for the 
values obtained with the forces NLl and NL3, respectively. 



58 



TABLES 



TABLE L The values of the parameters of different sets for the finite range Gogny forces. 



Parameter 


DIS (Ref. 




DIP (Ref. 


il) 


Dl (Ref. 11) 


fii (fm) 


0.7 




1.2 


0.9 




1.44 


0.7 


1.2 




-1720.30 




103.64 


-372.89 




34.62 


-402.4 


-21.30 


Bi 


1300.00 




-163.48 


62.69 




-14.08 


-100.0 


-11.77 


Hi 


-1813.53 




162.81 


-464.51 




70.95 


-496.2 


37.27 


Mi 


1397.60 




-223.93 


-31.49 




-20.96 


-23.56 


-68.81 
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TABLE II. The non-linear parameter sets NLl, NLSH and NL3. The masses are given in MeV, 
the parameter 52 in fm~^, while the rest of the parameters are dimensionless. The nuclear matter 
properties, predicted with these effectives forces, namely, the baryon density po (in units fm~^), the 
binding energy per particle E/A (in MeV), the incompressibility K (in MeV), the effective mass m*/m 
and the asymmetry parameter J (in MeV) are also shown. 



Parameter 


NLl 


NLSH 


NL3 


Masses 


rriN 


938.0 


939.0 


939.0 




492.25 


526.059 


508.194 




795.36 


783.0 


782.501 


rup 


763.0 


763.0 


763.0 


Coupling constants 


9a 


10.138 


10.4444 


10.217 


92 


-12.172 


-6.9099 


-10.431 


93 


-36.265 


-15.8337 


-28.885 


9oj 


13.285 


12.945 


12.868 


9p 


4.976 


4.383 


4.474 


Nuclear matter properties 


Po 


0.153 


0.146 


0.148 


E/A 


-16.488 


-16.346 


-16.299 


K 


211.29 


355.36 


271.76 


m* /m 


0.57 


0.60 


0.60 


J 


43.7 


36.1 


37.4 
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TABLE III. The spin assignment for the lowest state Iq of the unhnked or tentatively linked 
ji92p|-|j jj.g^i gj3 bands. The lowest transition energies E^{Iq + 2 ^ Iq) are also shown. 



Nucleus 



Ref. 



E^{Io + 2^Io) [keV] 







190jjg 
192Hg 
192pb 
196pb 
198pb 
198po 



m 
m 

\m 



316.9 
214.4 
214.8 
171.5 
305.1 
175.9 



12+ 
8+ 
8+ 
6+ 
12+ 
6+ 
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